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Abstract 



We calculate the heavy quark potential from the magnetic current due to monopoles in 
four dimensional SU(2) lattice gauge theory. The magnetic current is located in con- 
figurations generated in a conventional Wilson action simulation on a 16 4 lattice. The 
configurations are projected with high accuracy into the maximum abelian gauge. The 
magnetic current is then extracted and the monopole contribution to the potential is 
calculated. The resulting string tension is in excellent agreement with the SU(2) string 
tension obtained by conventional means from the configurations. Comparison is made 
with the U(l) case, with emphasis on the differing periodicity properties of SU(2) and 
lattice gauge theories. The properties of the maximum abelian gauge are discussed. 



I. INTRODUCTION 



In this paper, we report on our calculations of the SU(2) string tension using mono- 
poles. The monopoles were located in SU(2) lattice gauge theory configurations after 
the configurations were projected with high accuracy into the maximum abelian gauge. 
These lattice calculations were motivated by an approach to continuum confinement out- 
lined by 't Hooft some time ago In 't Hooft's framework, the first step is a partial 
gauge-fixing, applied only to those gauge fields which are "charged" , or have off-diagonal 
generators in the Lie algebra of the gauge group. The central idea is that the monopoles 
associated with the abelian gauge invariance left after partial gauge-fixing will control 
non-perturbative phenomena. This abelian gauge invariance is associated with the fields 
whose generators are diagonal, called "photons". For an SU(2) gauge group with the 
generator T 3 diagonal, the gauge field A 3 ^ is the abelian field or photon, and gauge-fixing 
is done only on A a „, a = 1, 2, or equivalently the charged fields 

For the particular choice of gauge-fixing condition known as the maximum abelian 
gauge, the continuum functional 

G c = \nl ((4) a + (4) 2 ) ^ = E / <wtw-)* x (i) 

is minimized over all SU(2) gauge transformations, leading to the conditions 

(d, + igAl)w; = (d, - igAl)W- = 0, 

where g is the SU(2) gauge coupling. In the remainder of this Section, we discuss how 
this is turned into a specific calculational scheme on the lattice. Our results are described 
in Section |I| . Section |TTTj contains a discussion of the maximum abelian gauge and our 
conclusions. 
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A. Lattice Gauge-Fixing 



The 577(2) lattice gauge theory is built out of link variables U^x), 



U^x) = e igaA ^, 



where f = a/2 are the generators of 577(2) in the fundamental representation, and a is 
the lattice spacing. On the lattice, the maximum abelian gauge is obtained by maximizing 
the lattice functional 



over all 577(2) gauge transformations ||. It is easy to show that in the continuum limit, 
maximizing Gi is equivalent to minimizing G c . 

At the maximum, Gi will be stationary under gauge transformations. The demand 
that Gi be stationary with respect to a gauge transformation at an arbitrary site y leads 
to the requirement that 



be diagonal. This can be accomplished by a gauge transformation However, the 

value of X at the nearest neighbors of y is affected by Q{y), so the diagonalization of X 
over the whole lattice must be done iteratively. 

After gauge-fixing, there is still manifest U(l) gauge invariance, and it is useful to 
factor a U(l) link variable from U^(x), writing U^(x) = u^{x)w^{x), where u^(x) = 
exp(i0^r 3 ), and w^ix) = exp(i9^ ■ r), with 9 Z = 0. The U(l) gauge transformation 
properties of and follow upon applying an abelian transformation Q 3 = exp(ia(x)r 3 ) 



G i = E^T Ul(x)a 3 U^(x)a 3 



(2) 



X{V) = E [u,(y)a 3 Ul(y) + Ul(y - fi)a 3 U{y - /*) 



(3) 



to Up. 



u 



, M (x) -> Ql(x + fi)Uf,(x)fl 3 (x) 



w^(x) -> Ql(x)w^(x)Q 3 (x), 
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so transforms as a charged chiral field at x. The angle 0^ can be extracted from the 
matrix elements of by expanding the gauge-fixed SU(2) link in Pauli matrices, 
writing 

U» = U° + zJ2K-a k . (4) 

k=i 

Then 0j = 2 ■ arctan(C/^/[/°). 

— # 

The maximum abelian gauge globally suppresses or equivalently, tries to force 
to the identity matrix. Even so, it is non-trivial to expect that long range effects 
associated with confinement are totally isolated in u^. The first concrete calculations to 
test this were performed by Suzuki and Yotsuyanagi 0. In Wilson loops, they replaced 
each SU(2) link variable by the U(l) link variable and found that full SU(2) results 
were obtained for Creutz ratios. This did not work for forms of partial gauge-fixing other 
than the maximum abelian gauge. 



B. Monopoles in U(l) 

In U(l) lattice gauge theory, the usual form of a Wilson loop involves a line integral 
of the U(l) link variable (^^(x) taken around a path specified by an integer-valued current 

(W u(1) ) = (exp(iJ2<t>lJ^y (5) 

where (■) denotes the expectation value over the ensemble of U(l) configurations. In 
addition, it is well established that confinement in U(l) is via monopoles A Wilson 
loop originally expressed as in Eq.([|) can be factored into a perturbative term arising 
from one photon exchange, times a non-perturbative term arising from monopoles,Q 



*Eq.(||) can be derived as an exact formula only for the Villain form of the U(l) action. However, 
in Ref. it was shown to work for other forms of the action, provided a coupling constant 
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(W U{1) )= (W phot ) ■ (W mm ) . (6) 

An explicit formula for (W mon ) in U(l) is obtained by writing J p as the curl of a Dirac 
sheet variable ||; = d u D^ u , where d u denotes a discrete derivative. Then (W mon ) is 
given by 



{W mon ) = 




where (•) denotes the sum over configurations of magnetic current. The sheet variable 
Dp V is not unique. For the usual case of an R x T loop with \J^\ = 1, a useful choice is 
to set D^ v = 1 on the plaquettes of the flat rectangle with boundary J M , and D^ v = on 
all other plaquettes. In Eq.(^), F* is the dual of the field strength due to the magnetic 
current; F*(x) = \e^ ua pF a p{x). The field strength itself is derived from a magnetic vector 
potential A™, F^ = d^A™ - d u A™, where 

A™(x) = J2v(x-y)m,(y), (8) 

y 

and m M is the integer- valued, conserved magnetic current. 

In Ref. [Q, results characteristic of confinement in SU(2) were obtained by replacing 
the full Wilson loop by a U(l) loop expressed in terms of the link variable <p z . Since 
confinement in U(l) lattice gauge theory itself is via monopoles, this raises the possibility 
of a monopole explanation of confinement in SU(2). 

C. Periodicity in SU(2) 

Before applying Eqs.(^) and (H) to SU(2), account must be taken of the difference 
in periodicity for the standard U(l) and SU(2) actions. The above discussion was for a 
£7(1) action periodic in link angles with period 2ir. The SU(2) action is periodic in the 



mapping was used to calculate {W p hot}- 
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£7(1) link angle 0^ with period 4tt. This follows from the formula u^{x) = exp^^T^) , 
where t% = 03/2. Alternatively, the process of temporarily setting to the identity on 
every link transforms the ST/ (2) action into a £7(1) action with period 47r: 

where <fp = - d^ . 



S — f £ (l-cos(^)), (9) 



While Eq.(^|) is not intended as an approximation to the full SU(2) action, it does 
correctly reveal the monopole charges which will occur in 577(2). A Dirac string occurs 
when the plaquette angle (jp is an integer multiple of Air, rather than 2tt, so the magnetic 
current m M here is an even integer, or in other words the abelian monopoles which occur 
in S£7 (2) are Schwinger monopoles J7|,Q . 

The replacement of by the identity in the action as in Eq. @ is never actually done. 
The configurations are generated using full SU(2) dynamics. However, after projecting 
the configurations into the maximum abelian gauge, the are set to the identity on every 
link in the calculation of Wilson loops. Only the £7(1) link variable u^{x) is retained, so 
that the 5£7(2) Wilson loop becomes a £7(1) loop: 

1 



W SU{ 2) — ► - x 



exp (i (pp^j + exp (-i ^ c/^J^j 



(10) 



where the conserved line current J M has | J M | = |. Since in this approximation, the Wilson 
loop is built out of £7(1) variables which involve only 0^, perturbative exchange of gluons 
coupled to T\ and r 2 has clearly been suppressed. By analogy with the situation in £7(1), 
perturbative exchange of the neutral gluon or "photon" coupled to t% is still allowed, 
but is expected to reside in the 5£7(2) analog of the (W p hot) factor of Eq.(|), whereas 
the confining part of the potential is expected to reside in the factor {W mon ). That the 
confining potential resides solely in (W mon ) is a postulate which will be justified by our 
results. 

Since |J M | = ~, the calculation of (W mon ) for 5£7(2) involves a Dirac sheet variable 
with \D^ U \ = ~, along with a magnetic current which is an even integer. Both of 
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these arose from the 47r periodicity of SU(2) in the link angle <ft \. It is straightforward to 
transform back to the familiar case of 2n periodicity. Define (p^ = <pf l /2 = arctan(C/^/C/°), 
and require ^ G (— ir, it]. The location of the magnetic current starts with plaquette 
angles <M„ constructed from 

4>IM = 9^1 - d^l = + + A) - 4>%x + *>) - <^(x) (11) 

The plaquette angle <j? is resolved into a Dirac string contribution, plus a fluctuating 
part: 

$% = 27rn tlv + 4>%, (12) 

where 0^ G (— 7r,7r] and is an integer f8|. The integer- valued magnetic current 
fh^ = m^/2 is determined by the net flux of Dirac strings into an elementary cube: 

% = ~^^uaf3 ■ d u n a p. (13) 

We now simply use m M in Eq.(||) to obtain A^, and from it, F*. Finally (W mon ) for 
SU(2) is obtained from Eq.(|7]) where the Dirac sheet variable now has \D^ V \ = 1. We 
have in effect used the Dirac condition |J on the product of electric and magnetic charge 
to transform from An to 2tt periodicity . In the interaction of a distribution of magnetic 
current with a single electric charge, if the magnetic current is halved and the electric 
current is doubled, the same result is obtained, since the interaction only depends on the 
quantized product of electric and magnetic charge. 

D. Monopole Location in 1-Cubes 

To summarize, we use Eq.(|7|) to calculate (W mon ) for SU(2), after extracting the 
magnetic current m M from 577(2) configurations projected into the maximum abelian 
gauge. In the identification of m^, elementary cubes were used, so a magnetic charge is 
located at the center of a spacial 1-cube, likewise for other components of the magnetic 
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current. This is done mainly for practical reasons; if (say) 2-cubes were used instead, the 
effective lattice size would become 8 4 instead of 16 4 . Although this practical reason was 
dominant in our calculations, it is worth pointing out that there is no conflict between 
monopoles being extended objects, and using 1-cubes to locate them. The monopole 
location procedure finds a monopole by locating the end of a Dirac string. Therefore, we 
must be in a gauge where there are Dirac strings attached to monopoles. Consider for 
the moment the case of a 't Hooft-Polyakov monopole in the continuum ||10|| . The gauge 
where a Dirac string is present is usually described as having the Higgs field along the 
isotopic 3-axis. A completely equivalent description is that the gauge fields are in the 
maximum abelian gauge . In this gauge, the photon field A 3 ^ takes the form appropriate 
for an elementary point monopole. The extended structure of the monopole involves only 
the charged gauge fields, [IT]. In a gauge with these properties, the center of the 
monopole may be located by finding the end of its Dirac string, using the finest spacial 
scale available. 

In calculating the effects of monopoles on heavy quarks, the assumption is made that 
despite their finite size, only the long range Coulombic fields produced by monopoles 
contribute to the confining potential. This produces a clear calculational procedure for 
computing (W mon ), which as will be seen in the next Section, works very well. However, 
further theoretical understanding of this assumption is certainly needed. 



II. MONOPOLE CALCULATIONS 



A. Simulation and Gauge-Fixing 



Our simulations were done on a 16 4 lattice, using the standard Wilson form of the 
SU(2) action. Three (3 values were used; (3 = 2.40,2.45, and 2.50. At each /3, after 
equilibration, 500 configurations were saved. Saved configurations were separated by 20 
updates of the lattice, where a lattice update consisted of one heatbath sweep [0, plus 
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one or two overrelaxation sweeps [13]. Each of these configurations was then projected 
into the maximum abelian gauge using the overrelaxation method of Mandula and Ogilvie, 
with their parameter uj = 1.70 |TJj]. The overrelaxation process was stopped after the off- 
diagonal elements of X(x) of Eq.(Q) were sufficiently small. Expanding X(x) in Pauli 
matrices, 

3 

X = X\x)+iY,X\x)-a kl (14) 

k=l 

we used 

(l**| a )^E(l*W + l*W) (15) 

X 

as a measure of the average size of the off-diagonal matrix elements of X over the lat- 
tice, and required (jX ch \ 2 ^ < 1CT 10 . This condition was reached in approximately 1000 
overrelaxation sweeps. 

From each gauge-fixed SU(2) link, the U(l) link angle <^ was extracted using the 
formula ^ = arctan(L^/C/°), as described in Section |. Then making use of the plaquette 
angles (f> 3 , the magnetic current m M (x) was found. Although the values m M = ± 1 ± 2 
are allowed by Eqs.(|T2"D and (|TBJ), the overwhelming fraction of links carrying current had 
|m M | = 1; [m^l = 2 rarely occured. Only a few percent of the links actually carried 
current. We define f m as the number of links with non-zero current, divided by the total 
number of links, 4 • L 4 . Our results for f m are recorded in Table I, and agree well with 



those in the literature 15.16 



B. The Quark Potential from Monopoles 

With the magnetic current in hand, the monopole Wilson loops for 577(2) are 
completely determined. Applying Eqs.(|7|) and (|8|), for (3 = 2.40 and 2.50 all R x T loops 
up to 7 x 10 were measured and averaged over configurations. For (3 = 2.45, the maximum 
size was increased to 8 x 12. From the monopole Wilson loops, monopole contributions 
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to the potential, denoted by V mon (R), were extracted by performing straight line fits to 
ln((W(R,T) mon )) vs T. The fits for R > 2 were over the interval T — R+ 1 to T max , 
except for R max , where T = R to T max was used. Finally, the monopole contribution to the 
string tension a was extracted by fitting V mon (R) to the form V mon (R) = a/R + a ■ R + Vo, 
over the interval R = 2 to R = R m ax- The results are shown in Fig. (1), and tabulated in 
Table II. 

As a glance at either Table II or Fig. (1) shows, the monopole potentials are essentially 
linear at all values of R, with negligible Coulomb terms. This is expected on the basis of 
the discussion of Section |. The steps of first setting = 1 on each link of the Wilson 
loop, and then calculating only the monopole factor (W mon ) of the resulting U{1) loop, 
effectively suppresses the Coulomb terms arising from single gluon exchange. 

The crucial question for the present paper is whether or not the monopole string 
tensions agree with those for full SU(2). Of the various ways of determining the SU(2) 
string tension, the "torelon" method of Michael and Teper is perhaps the best comparison. 
The torelon method involves the temporal correlation between loops which are wrapped 
around the entire lattice in the spacial direction. The method determines the string 
tension directly, without the need to separate linear from Coulomb terms. For yfa, Michael 
and Teper give 0.258(2) at /3 = 2.4 on a 16 4 lattice and 0.185(2) at /3 = 2.5 on a 20 4 
lattice jnj. Both numbers are in excellent agreement with the monopole results of Table 
II. As an additional check that full SU(2) results were being obtained, we performed 
linear-plus- Coulomb fits to our own SU(2) Wilson loops, which were obtained from the 



same set of configurations using an analytic form of the multi-hit method [18|, along with 



the smearing method fll9[ , as noise reduction techniques. The resulting string tension is 
presented in Table III, and again there is agreement with the monopole string tension 
to within statistical errors. Statistical errors were estimated using standard jackknife 
methods. Similar results for the SU(2) string tension using monopoles have recently been 
obtained by Shiba and Suzuki [EO] . The Coulomb terms coming from the monopoles and 
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full SU (2) naturally differ, as explained above. No attempt to calculate the SU (2) analog 
of the factor {W p hot) is made in this paper. 



C. The String Tension and Magnetic Current Loop Size 

To further investigate confinement via monopoles, the magnetic current was resolved 
into individual loops, each of which separately conserves current. Not all sizes of loops 
are expected to contribute to the confining potential. For example in U(l) lattice gauge 
theory, it is possible to show analytically that a random distribution of loops much smaller 
in size than the Wilson loop under consideration does not affect the long range potential, 
but only renormalizes the Coulomb term. At the opposite extreme, again in U(l), it is 
known that the string tension can be calculated accurately using only the contribution 
coming from very large loops of magnetic current ET] . 



For the present case of SU(2), we first obtained a rough measure of how the current is 
distributed over loops of various size. This was done by counting the number of current- 
carrying links residing in loops of size up to and including 10, 20, 50, and 100 links. In 
Table IV, the results are shown as a fraction of the total current. As can be seen from 
the last column of the table, including all loops of magnetic current with up to 100 
links accounts for approximately half of the current-carrying links at each coupling. The 
remainder of the current consists of a small number (~ 3) of loops, each containing 
typically several hundred links. 

To study how the string tension depends on magnetic current loop size, we computed 
(W mon ) with a cut on loop size. Since experience with U(l) shows that large loops are 
what is important in confinement, in computing (W mon ) from Eq.(|7]), only those loops 
of magnetic current with more than n cut links were included. The results for the string 
tension with n cut = 50, 100, and 200 are shown in Table V. For n cut = 50 the answers 
are within statistical errors of the full string tension for all couplings. However, as n cut is 
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increased the string tension decreases steadily so that by n cut = 200, there is a statistically 
significant deviation from the full string tension. This shows that in SU(2) lattice gauge 
theory, the string tension cannot be explained by retaining only the very largest loops 
of current, unlike the situation in U(l) . To make this point quantitatively, we compare 
to our work in £7(1) lattice gauge theory on a 24 4 lattice. In Ref. [^lJI , we considered a 
£7(1) coupling corresponding to a string tension of a = 0.058(2), intermediate between the 
SU(2) string tension at (3 = 2.40 and f3 = 2.45. For this case and other nearby couplings, 
the £7(1) string tension is stable under an increase of n cu t to at least n cut = 1, 000, whereas 
deviations are already significant in 577(2) for n cut = 200. Admittedly, two different lattice 
sizes are being compared here, but since our 16 4 577(2) string tension is within statistical 
errors of 24 4 577(2) numbers, it is reasonable to assume that the distribution of loops we 
found on 16 4 would be similar to that on a 24 4 lattice. 

Small loops of magnetic current play no role in confinement for either 577(2) or £7(1). 
In the present work, we performed additional calculations of (W mon ) with a cut on loop 
size, but this time including only loops of current with less than n cut links. The string 
tension was statistically zero for n cut = 50 and 100 for all three values of (3. The Coulomb 
terms were small, of magnitude ~ 10% of the Coulomb term for full SU (2), and attractive 
in sign. For n cut = 200, a very small contribution to the string tension was seen for (3 = 2.5. 

We may summarize our investigation of current loops by saying that for SU(2), some- 
what more than half the magnetic current plays a role in the confining potential. The 
string tension can be explained by retaining all magnetic current loops of intermediate 
size (~ 50 links) and larger. An reasonable guess for how large a loop must be to play a 
role in the confining potential would be a physical extent 0(1/ ^/a) or greater. For U(l), 
there is a class of huge magnetic current loops present only in the confined phase, and 
the string tension can be explained with a restriction only to these very large loops. 
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D. Abelian Wilson Loops 



Another possible route to the string tension is Eq. (|10D , which expresses the assumption 
that after partial gauge-fixing, the physics of confinement is contained in the £7(1) link 
variable <fP. So far, we have concentrated on the extraction of the string tension from 
(W mon ), which is based on the further assumption that the £7(1) approximation to the 
full £77(2) Wilson loop, Eq . (]T0|) , can be factored into a photon part (W p hot), times a 
monopole part (W mon ). The physical motivation for this was to investigate the monopole 
confinement mechanism, but the method also has clear computational advantages. The 
noise level in the monopole Wilson loops is much lower than in the full £77(2) loops, 
and the string tension is easier to identify, due to the very small Coulomb terms in the 
monopole contribution to the potential. For completeness, we have returned to Eq. (p70|) to 
calculate the £7(1) approximation to the ££7(2) Wilson loops directly in terms of the link 
variable 0^, and extracted the string tension and Coulomb term. The statistical accuracy 
of these results is poor compared to those obtained from the full SU(2) Wilson loops, due 
to the absence of the multi-hit technique. Partial gauge- fixing produces a sequence of £7 (1) 
configurations, but the £7(1) action that would produce these configurations is unknown. 
Lacking the usual noise reduction techniques, we perform a simple average over these £7(1) 
configurations to obtain Wilson loops, and extract the string tension and Coulomb term 
from the resulting potential. The Coulomb term is now significant, but still less than that 
for full ££7(2). This is expected since exchange of neutral gluons or photons is allowed 
in Eq . (|10|) . The string tension is consistent with the full ££7(2) string tension as well as 
that extracted from {W mon ) |p2]| . In Fig. (2), we show the potentials derived from the 
full ££7(2) loops, from the monopole loops, and from Eq.(^), for (3 = 2.45. A constant 
has been added to the monopole and £7(1) potentials for the purpose of comparing the R 
dependence of the three potentials in Fig. (2). The shift of the potentials by a constant 
is not physically relevant in that it does not affect the functional dependence on R of the 
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potentials. The plots for (3 = 2.40 and 2.50 are similar. 



III. CONCLUSIONS AND SUMMARY 

The present work was done entirely in the maximum abelian gauge, stimulated by 
the pioneering work of Suzuki and Yotsuyanagi ||. In an average way over the lattice, 
the maximum abelian gauge forces the fluctuations in the charged sector to be as small 
as possible. This gauge also has a variational formulation, (Eqs.(|l]) and (|^)), and is 



renormalizable in the continuum limit |23]. While these are desirable features, it should 
be possible to capture the physics of confinement with other forms of partial gauge-fixing. 
That this has not been possible so far in lattice calculations may have to do with the 
monopole location procedure. As mentioned in Section [1D| , the location procedure finds 
a monopole by locating the end of its Dirac string, so if monopoles are to be correctly 
located using the abelian flux over 1-cubes, the gauge which is used must attach a string 
to a monopole at distance scales right down to the lattice spacing. For a 't Hooft-Polyakov 
monopole in the continuum, the requirement that the string go all the way to the center 
of the monopole singles out the maximum abelian gauge. It is easy to check that another 
gauge discussed by 't Hooft [p]], where the charged field strength Gf 2 is set to zero, does 
not have this property. Starting from the "hedgehog" form of the monopole solution, 
going to the gauge where Gf 2 = will be accomplished by the same gauge transformation 
as going to the maximum abelian gauge in the region far from the monopole, so here 
there is no distinction between the two gauges. However, inside the extended structure 
of the monopole, the two begin to deviate and only the maximum abelian gauge has a 
Dirac string extending to the origin of the monopole.! Assuming a similar phenomenon 
occurs on the lattice, a calculation which uses 1-cubes to locate the monopoles in the 



^At the center of the monopole, the hedgehog solution already satisfies Gf 2 = 0. 
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lattice version of the gauge Gf 2 = may then locate the monopoles incorrectly. This 
interpretation of the difference between the maximum abelian gauge and other forms 
of partial gauge-fixing is supported by work on three-dimensional SU(2) lattice gauge 
theory by Trottier, et al. p4[ . They show that the density of monopoles as found in 
the maximum abelian gauge has the correct scaling law as a function of coupling, using 
1-cubes to locate the monopoles. However, for other forms of partial gauge-fixing such 
as Gf 2 — 0, the correct scaling law is obtained only after the size of the cube used to 
locate the monopoles is increased considerably. This would say that while other forms 
of partial gauge-fixing are in principle equivalent, for practical reasons having to do with 
lattice size the maximum abelian gauge is likely to continue to be favored in calculations 
with monopoles. 

We have presented evidence in favor of 't Hooft's picture of confinement. Although 
our lattice size and number of configurations is moderate, the numbers obtained from 
monopoles for the string tension are in excellent agreement with full SU (2) numbers on 
larger lattices with higher statistics. It would be of interest to repeat the present work on 
a larger lattice with several thousand configurations, in order to search for any possible 
systematic difference between the full SU(2) string tension and the string tension deduced 
from monopoles. A larger lattice would also allow a move toward weaker coupling and a 
smaller string tension. As the correlation length defined by the string tension grows, so 
presumably must the physical size of monopoles which has been ignored in the present 
work. An interesting issue is whether the string tension will continue to be dominated by 
the long range Coulombic fields of the monopoles. Study of the gluon propagator offers 
a different line of attack on this question. We are presently calculating both charged and 
neutral gluon propagators from our 16 4 configurations, after a final gauge-fixing which 
puts the photon link variable 0^ in the lattice Landau gauge [p5|| . A detailed report 
will be presented elsewhere, but our preliminary calculations show that the neutral gluon 
or photon propagator is intrinsically larger and of longer range than the charged gluon 
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propagator. 
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CAPTIONS 



Table I. The fraction of links carrying magnetic current at each value of (5. 

Table II. The string tension and Coulomb coefficient obtained from fits to the potentials 
calculated using the monopole Wilson loops at each (3. 

Table III. The string tension and Coulomb coefficient for the full SU(2) potential at each 
P- 

Table IV. The fraction of total magnetic current contained in loops with < 10,20,50, and 
100 links. 

Table V The monopole string tension from all current loops with > 50,100, and 200 
links. 

Figure 1. The potentials extracted from monopole Wilson loops at f3 — 2.40 (circles), 
2.45 (squares), and 2.50 (triangles). The solid lines are the linear-plus-Coulomb fits to 
each potential. 

Figure 2. Comparison of the monopole potentials (circles), the U(l) approximation to 
full SU(2) potentials (squares) and the full 577(2) potentials (triangles) at f3 — 2.45. 
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ft 

p 


r m 


2.40 


2.75(6) x 10~ 2 


2.45 


1.95(1) x lO- 2 


2.50 


1.36(1) x 10~ 2 



Table I 



p 


a 


a 


2.40 


0.068(2) 


0.01(1) 


2.45 


0.051(1) 


0.02(1) 


2.50 


0.034(1) 


0.01(1) 



Table II 



p 


a 


a 


2.40 


0.072(3) 


-0.28(2) 


2.45 


0.049(1) 


-0.29(1) 


2.50 


0.033(2) 


-0.29(1) 



Table III 
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max loop size 


p 


10 


20 


50 


100 


2.40 


0.33 


0.41 


0.46 


0.48 


2.45 


0.38 


0.46 


0.51 


0.54 


2.50 


0.43 


0.51 


0.56 


0.59 



Table IV 



p 


<j(> 50) 


cr{> 100) 


cr(> 200) 


2.40 


0.068(3) 


0.066(3) 


0.063(3) 


2.45 


0.050(3) 


0.049(3) 


0.044(3) 


2.50 


0.035(2) 


0.034(2) 


0.029(2) 



Table V 
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Figure 1 
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Figure 2 
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